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We analyze the optical properties of one-dimensional (ID) "PT-symmetric structures of arbitrary 
complexity. These structures violate normal unitarity (photon flux conservation) but are shown to 
satisfy generalized unitarity relations, which relate the elements of the scattering matrix and lead to a 
conservation relation in terms of the transmittance and (left and right) reflectances. One implication 
of this relation is that there exist anisotropic transmission resonances in "PT-symmetric systems, 
frequencies at which there is unit transmission and zero reflection, but only for waves incident from 
a single side. The spatial profile of these transmission resonances is symmetric, and they can occur 
even at PT-symmetry breaking points. The general conservation relations can be utilized as an 
experimental signature of the presence of PT-symmetry and of "PT-symmetry breaking transitions. 
The uniqueness of PT-symmetry breaking transitions of the scattering matrix is briefly discussed 
by comparing to the corresponding non-Hermitian Hamiltonians. 
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I. INTRODUCTION 

Motivated by fundamental studies of T-'T-symmetric 
quantum Hamiltonians [IHS], PT-symmetric photonic 
structures have attracted considerable interest in the past 
few years. These are structures with balanced gain and 
loss; in the case of a one-dimensional (ID) structure, this 
means that there is a symmetry point (chosen to be the 
origin, a; = 0) around which the linear index of refraction 
satisfies n*{—x) = n{x). Such structures were first stud- 
ied in Refs. [H Ej and were shown to exhibit a variety of 
exotic photon transport phenomena, such as double re- 
fraction [S], power oscillations [SH7], and non-monotonic 
behavior of the transmission loss with increased dissipa- 
tion [8^ . The initial studies focused on parallel waveguide 
structures with alternating loss and gain, in which the 
transverse variation of the electrical field, in the parax- 
ial approximation to the wave equation, maps precisely 
onto a ID or discrete Schrodinger equation, similar to 
the earlier quantum studies [IHS] . The parallel waveguide 
realization of VT symmetric photonic structures has re- 
cently found a promising application to compact optical 
isolators and circulators |T0| . 

Recently, several authors have studied PT-symmetric 
cavities and heterostructures [HHS], as well as general 
VT scattering systems jT3], using the full scalar wave 
equation, in the case that it obeys at least one VT sym- 
metry operation. The current authors in particular em- 
phasized the existence in such systems of PT-symmetric 
and PT-broken phases of the electromagnetic scattering 
matrix (S-matrix). For the ID case, the eigenvalues of 
the S-matrix are unimodular in the PT-symmetric phase, 
as they are in unitary systems, but photon flux is not 
conserved for most scattering processes, whereas in the 
PT broken phase the S-matrix eigenvalues have recip- 



rocal magnitudes, one greater than unity (corresponding 
to ampliflcation) , and the other less than unity (corre- 
sponding to attenuation) . We and others pTHT5] pointed 
out the existence of novel singular points in the broken 
symmetry phase, which we refer to as CPA-laser points. 
At these points one of the S-matrix eigenvalues goes to 
infinity (the usual lasing threshold condition), while the 
other goes to zero. The latter phenomenon corresponds 
to coherent perfect absorption (CPA)[Tni |T7]: in which 
a specific mode of the electromagnetic field, the time- 
reversal of the lasing mode, is completely absorbed. For 
PT-symmetric structures, these two phenomena must co- 
incide [mils]; i-e. at the laser threshold, in addition to a 
radiating mode of self-oscillation, there always exists an 
incident field pattern which, instead of being amplified, 
is completely attenuated. 

The rich behavior of ID PT-symmetric photonic struc- 
tures violates the standard intuition that optical struc- 
tures can be characterized by their single-pass gain or 
loss, which is always zero in these systems. The coinci- 
dence of both lasing and perfect absorption, and more 
generally the reciprocal amplification and attenuation 
displayed by the S-matrix eigenvalues, is a strict con- 
sequence of the symmetry property of the S-matrix for 
such structures. In Ref. [13j this was expressed in arbi- 
trary dimensions by the relation 

{VT)Siu*){VT)^S-'{ij), (1) 

where P is the parity operator (or indeed any discrete 
symmetry operator with P^ = 1), and T is the time- 
reversal operator (in the representation we will employ, 
this can be taken as the complex conjugation operator). 
By comparison, a T-symmetric unitary S-matrix would 
obey TS{u;*)T = S~\u;). 

The set of S-matrices obeying Eq. (fTl) can be shown to 
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be isomorphic to a pseudo-unitary group, which in the 
ID case is just f/(l, 1) [18]. In physical dimensions higher 
than one, there can be more than two input and output 
channels, and it is possible for the S-matrix to be in a 
mixed "phase" with one subset of the eigenvalues form- 
ing "T'T-broken" amplifying/attenuating pairs and the 
remaining eigenvalues being "PT-symmetric" and flux 
conserving. For ID structures, however, there are only 
two eigenvalues, and they must either be both unimodu- 
lar, or a non-unimodular inverse conjugate pair — except 
at the 'PT-transition point, an exceptional point at which 
the S-matrix has only one eigenvector and eigenvalue |13j . 

Several specific cases of ID PT-symmetric structures 
have been studied [TTUISl fT9] and in addition to the in- 
teresting CPA-laser behavior, other intriguing proper- 
ties have been found, such as unidirectional invisibility 
[19j . It is thus worthwhile to see what specific proper- 
ties T^T-symmetry imposes on transmission and reflec- 
tion in arbitrary VT structures, in both the symmetric 
and broken-symmetry phases. That is the goal of the 
current work. In Section |llj we show that ID VT struc- 
tures obey certain strong conservation relations, which 
could be employed experimentally to determine if a given 
structure has realized VT symmetry. In Section III we 



examine a consequence of these conservation relations: 
the existence of transmission resonances in which the re- 
flectance vanishes only for waves incident from one side of 
the structure, which we refer to as anisotropic transmis- 
sion resonances (ATRs). The uni-directional invisibility 
phenomenon found by Lin et al. [19| is a special case of 
these ATRs. In Section [TV] we derive a separate relation 
for the boundary between the PT-symmetric and VT- 
broken phases of the S-matrix, involving the reflectance 
and transmittance for one-sided scattering processes. In 
Section [Vj we show that our conventional definition of 
the S-matrix and its eigenvalues is physically meaning- 
ful, and in particular that its phase boundary can be 
related to PT-breaking transitions in the spectrum of 
some T'T-symmetric Hamiltonian. 



II. GENERALIZED UNITARITY RELATIONS 

We begin, following Longhi [12], with the ID transfer 
matrix M, defined by (see Fig. [T]): 



M 



(2) 



For a T^T-symmetric heterostructure, the components of 
M obey the following properties |12| : 

MM = M*,{u;*), Afi2(2i)(^) = -Mr2(2i)(^*), (3) 

where uj is the frequency of the incident/scattered beams. 
For real w, these relations imply M22 = Mi-^ and 
Re[Mi2] = Re[M2i] = 0, which enables us to parame- 
terize M as 



M = 



a* ib 
—ic a 



(4) 
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FIG. 1; Schematic of a ID "PT-symmetric photonic het- 
erostructure, consisting of an arbitrary number of layers that 
are PT-symmetric about a:: = 0, i.e. n{x) = n{—x)*. G (red) 
and L (gray) indicate gain and loss regions, and different color 
tones indicate different amplification/absorption strengths. 



It is determined by three independent real quantities, i.e. 
b and the phase and amplitude of a. The parameter c is 
related to |a|, 5 by 



be 



(5) 



which arises from the quite general condition det(M) = 1 
[2U] . The parametrization using a, 6 is valid except when 
-^12 = 0; in that case, \a\ — 1 and c replaces b as the 
third independent parameter. 

In the following discussion we assume non-vanishing 
Mil and M22, which holds everywhere except at CPA- 
laser points [H]. The S-matrix is defined by 



= S 



(6) 



where /r are the refiection coefficients for light incident 
from the left and right respectively, while t is the trans- 
mission coefficient, which is independent of the direction 
of incidence. The parametrization Q gives 



S 



ib 1 
1 ic 



(7) 



Thus the refiection coefficients are r^, = ib/a and = 
ic/ a, which are unequal in magnitude but can differ in 
phase by only or tt; and the transmission coefficient is 
t = 1/a. Note that S satisfies the symmetry relation ([I]), 
with V = (? 0) and T the complex conjugation opera- 
tor. Using (pi), we obtain the following exact "generalized 
unitarity relation" : 



ThrR 



(8) 



This leads to the conservation relation 



IT - 1 



- ^RlRr, (9) 

where Rl/r = ^L/flP are the two refiectances and 
T = l^p is the transmittance. In addition, Eqs. (TjSl lead 
to phase relationships among the refiection and transmis- 
sion coefficients: 



bR^^L, ifr<i; 

bR = ^L+T^, ifr>l; 
bL,R = 0t ± 7r/2, 



(10) 
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where (I)l,r, 4>t are the phases of the reflection and trans- 
mission coefficients. 

Eqs. ( |8|9p are the central results of this work. They 
are valid for all ID photonic heterostructurcs with VT- 
symmetry; two examples are shown in Fig. [2j 




FIG. 2; (a) Reflectance and transmittance of a ID VT- 
symmetric structure of index n = 2 ± 0.2i and length L. 
Rl, Rr,T are indicated by the red, green, blue curves. Zeros 
of the reflectances and the corresponding anisotropic trans- 
mission resonances [T = 1) are marlced by dotted grey lines. 
The quantity RlRr + 2T - T'^ = 1 is plotted as the dashed 
black line to demonstrate the conservation relation ([9|. (b) 
Phases of tl, tr,, and t in (a), demonstrating the reflection 
phase jumps at each ATR. (c) Same quantities are plotted as 
in (a), but the structure now has a passive region of length 
2L/5 in the center. For this structure we see that it is pos- 
sible to have "accidental" flux-conserving points at which 
Rl ~ Rr{= R ^ 0) and hence T + R — 1. Similarly, there are 
accidental pseudo-unitary points for which T ~ R = 1. They 
are indicated by vertical dotted and dashed lines, respectively. 



For T < 1, Eq. (|9|) becomes T + ^/RlRr = 1. This 
is an intriguing generalization of the more familiar con- 
servation relation R + T = 1, which applies to uni- 
tary (T-symmetric) S-matrices for which the left and 
right reflectances are necessarily equal. In the VT- 
symmetric case, the geometric mean of the two re- 
flectances, ^/RlRji, replaces the single reflectance R. 
Therefore, when T < 1, the scattering of a single inci- 
dent wave from one side of the structure is sub-unitary 
(some flux is lost) and the scattering from the other side 
is super-unitary (some flux is gained). As an exception, 
there can be an accidental degeneracy at which R^ = Rr, 



in which case the scattering from both sides conserves 
flux. Such special cases do occur as a continuous param- 
eter such as frequency is varied for non-trival VT systems 
(Im[n(x)] ^ 0), as shown in Fig. [2](c). 

For T > 1, all single-sided scattering processes are 
super-unitary, and the conservation relation (|9| can be 
re-written as T — ^/RlRh = 1. Accidental reflectance 
degeneracies {Rl — Rr) are also possible in this regime, 
giving the usual pseudo-unitary conservation relation 
T - R = 1, as shown in Fig. [^c). AU of these quan- 
tities actually diverge when approaching the CPA-laser 
points, but they still satisfying the conservation relation 

Finally, we see that for T — 1, one of the reflectances 
must vanish (the other typically does not). Hence, the 
scattering for that direction of incidence is flux conserv- 
ing, similar to resonant transmission in unitary struc- 
tures. This phenomenon is analyzed in greater detail in 
Section Uni 

Interestingly, the S-matrix describing three-wave mix- 
ing in the undepleted pump approximation corresponds 
to the the special case where Rl — Rr [221 123] . The case 
T + R = 1 describes frequency conversion by absorption 
or emission of a pump photon, and T — R = 1 describes 
parametric amplification of both signal and idler by down 
conversion of pump photons. The relevance of a special 
case of VT symmetry to optical parametric amplification 
and conversion has only very recently been appreciated. 

An experimental concern in all VT systems is how to 
confirm that one has truly realized a structure with VT- 
symmetry, i.e. that the gain and loss are balanced and 
the real index is symmetric. Eqs. (8p0 ) are strong con- 
straints on the allowed scattering processes with a single 
incident beam for VT systems, and can be used to test 
how close one is to the ideal symmetric structure. 



III. ANISOTROPIC FLUX-CONSERVING 
TRANSMISSION RESONANCES 

As we have noted, Eq. (9]) implies an interesting phe- 
nomenon: there exists a flux-conserving scattering pro- 
cess for incident waves on a single side if and only if 
T = 1, and one of Rl or Rn vanishes. We refer to such a 
process as an anisotropic transmission resonance (ATR). 
ATRs are different from the accidental flux-conserving 
processes that can occur for T < 1; those, as we have 
seen, are accessible from either direction of incidence 
{Rl = Rr)- ATRs are a generalization of the flux- 
conserving transmission resonances of unitary systems, 
which are independent of the incidence direction. In 
Fig. [3) we show how two ATRs evolve out of a single 
transmission resonance of the unitary system as balanced 
gain and loss is added. Within the same structure, ATRs 
can occur for both left and right incidence, as the fre- 
quency is varied, but generally at different frequencies 
(to occur at the same frequency, a "doubly accidental" 
degeneracy Rl = Rr = would have to occur, requiring 
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a second tuning parameter). 

A surprising property of ATRs is that their intensity 
profile is spatiaUy symmetric. This can be shown from 
the following analysis. If E{x) is the spatial profile of a 
lcft-/right-going transmission resonance, then by a VT 
operation E*{~x) is also a left-/right-going transmission 
resonance of the same structure. Since these two states 
happen at the same frequency, they must be identical (up 
to a phase 0) by the requirement of uniqueness: 



E*{-x)=e"''E{x). 



(11) 



Hence, the intensity satisfies I{x) = \E{x)\'^ = /(— x). 
This result is consistent with the intuitive expectation 
that in order to conserve flux the photons must on av- 
erage spend equal amount of time in the loss and gain 
regions of the structure. Except at the ATRs, intensities 
do exhibit asymmetry for single-sided incidence, and in 
particular this is the case for a wave incident from the 
side with non-vanishing reflectance; see Fig. [3|^e),(f). 

Fig. |3] shows two ATRs of a multi-layer structure, 
one for each incidence direction, occurring at differ- 
ent frequencies. The frequencies are very similar be- 
cause Im[n] is not very large and both ATRs arise 
from a bi-directional transmission resonance of the uni- 
tary (Im[n] — 0) heterostructure. As we add gain and 
loss to the unitary heterostructure, while preserving the 
■pT-symmetry, the transmission resonances separate and 
their spatial profiles become more distinct. Fig. 
shows the asymmetric intensity profiles for waves inci- 
dent at the ATR frequency, but from the "wrong" side 
(the side with non- vanishing reflectance). The asymme- 
try increases as the two ATRs move further apart with 
increasing gain/loss, as shown in Fig.jsj^f). 

Let us refer to the left and right halves of a VT- 
symmetric heterostructure as U, V. We can write the 
reflection and transmission coefficients coefficients for the 
whole structure (r^, rji, and t) in terms of the refiection 
and transmission coefficients for the U and V segments: 



rL = 



rL,u 



' L,U 



1 _ pi{au+av)r* r* 
^ ^ ' L,U' R,V 



rR.v 



t 



^ ^ ' L,U' R,V 



(12) 

(13) 
(14) 



Here, a^z/y = 2Arg[t[//v]- Note that if either vl^u = or 
rji v = at some w, corresponding to a transmission res- 
onance of U/V in right/left direction, the transmittance 
for the full structure will also be unity. 

Thus one type of ATRs can arise from resonances of 
either halfoi the VT system. This follows from VT sym- 
metry. First, using the time-reversal operation, a trans- 
mission resonance of S(nk) from the left must be a trans- 
mission resonance of S(n*k) from the right (interchange 
gain and loss regions and interchange incoming and out- 
going amplitudes) [17 . Second, the S-matrix of the right 
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FIG. 3: Asymmetric transmission resonances in a multi-layer 
heterostructure. (a) Transmittance and reflectances as a func- 
tion of frequency for Im[n] 7^ and, (b), for Im[n] = 0. 
Mode is a bi-directional transmission resonance in the lat- 
ter case. The structure has a constant Re[n] = 2 and con- 
sists of 50 layers with Im[n] increasing (decreasing) stepwise 
from 0.004 (-0.004) to 0.1 (-0.1) towards the center in the 
loss (gain) half. The frequencies are LtjL = 25.275, 25.139 in 
Mode 1 (wave incident from gain side) and 2 (wave incident 
from the loss side). Spatial profiles shown are: (c) the ATRs, 
(d) unitary bi-directional transmission resonance, (e) "wrong" 
side scattering at the ATR frequencies, and (f) "wrong" side 
scattering at Rr = in the same structure but with larger 
Max[Im[n]] = 0.17, showing a stronger asymmetry. 
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FIG. 4: Graphical explanation of a trivial ATR arising from 
the transmission resonance of the left half in the right direc- 
tion. VT symmetry requires that it also be a resonance of the 
right half for left incidence (see argument in text). Depending 
on the particular VT structure there may or may not be triv- 
ial ATRs; for example the simple heterostructure of Fig. [2]^a) 
has no such ATRs. The ATRs of primary interest are those 
that arise from multiple scattering between the left and right 
halves of the structure. 



hand side of a VT structure is VS{n*k), so the right half 
of the VT structure must have a resonance for waves 
incident from the left side as well, if its left side does. 
Therefore the composite structure will have an ATR if 
either half does (r^.;/ = or rn v = 0)- This argument 
is illustrated graphically in Fig. |4j we refer to these as 
trivial ATRs. 

ATRs also occur when Arg[rL,;7] or Arg[rij^y] equals 
{au + av)/^ and involve multiple scattering between the 
sub-units. It is straightforward to check that at such 
points T = 1 and Rl{Rr) = 0. It can be shown that a 
single layer of gain or loss in a lossless environment (e.g. 
in air) does not have transmission resonances in general, 
and we show in the appendix that all the ATRs in Fig. [2] 
are of this type and are thus "non-trivial" . 
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As already noted, for an ATR to be bi-directional, a 
doubly accidental degeneracy is needed either in the am- 
plitude of (7 and rny {^l.u — fR.v = 0) or their 
phase {Aig[rL,u] = Arg[r_R,y] = {au + av)/2). This is 
highly unlikely, unless one can tune an additional con- 
tinuous parameter other than the frequency, so in the 
generic case all transmission resonances of VT structures 
are uni-directional. 

In a recent work, Lin et al. [19 have studied a ID 
T'T-symmetric Bragg structure of alternating dielectric 
layers with appropriate gain and loss, and discovered an 
ATR with an additional property which they refer to 
as "uni-directional invisibility" . Not only do they find 
T = l,_Ri = 0, ^ (or vice versa), as dictated by 
Eq. they also find that at this ATR the transmis- 
sion phase (j)t — Q, corresponding to zero phase delay of 
the signal compared to free propagation. For this reason 
there would be no signature of the presence of the struc- 
ture in either amplitude or phase of the received wave, if 
the wave is sent from the correct side (there would be a 
signature of course in the reflected wave if sent from the 
wrong side). This condition, that = at the ATR, is 
not required by our generalized unitarity relations and is 
specific to their structure [21]. 

The existence of non-trivial ATRs is independent of 
whether the S-matrix is in the PT-symmetric or VT- 
broken phase |13j ; They can even occur at the symmetry- 
breaking exceptional point (see the following section). 
However, we do find for the simple gain-loss heterostruc- 
ture of Fig. [2f^a),(b) that the ATRs disappear soon after 
the lasing threshold is passed in the broken symmetry 
phase, since in the large Ini[n]fcL limit R^Ri^ approaches 
unity asymptotically. This is not the case for more com- 
plicated VT structures such as that of Fig. . The dif- 
ferent behaviors of the two cases are illustrated in Fig. [TO] 
of the Appendix and its origin is discussed. 



IV. PHASE TRANSITION BOUNDARIES 

A ID T'T-heterostructure can undergo a spontaneous 
symmetry-breaking transition in the eigenvalues and 
eigenvectors of its S-matrix, as either w is increased at 
fixed gain/loss or as gain/loss is increased at fixed a; [T5]. 
In the symmetric phase, the VT operation maps each 
scattering eigenstate back to itself, whereas in the broken 
symmetry phase each scattering eigenstate is mapped to 
the other. At the symmetry breaking exceptional point, 
there is only one eigenvector and so both cases coincide. 

Let Ai_2 be the eigenvalues the S-matrix of a VT- 
symmetric heterostructure and Vi 2 be the ratios of the 
two amplitudes of the corresponding eigenstates. It fol- 
lows from the S-matrix parameterization ([t]) that 



I 

2a 

i 

2 



(6 + c) ± v/(6-c) 
(c-6)±y(6^ 



2-4 



(15) 
(16) 



These equations imply that A1A2 = J^ii^2 = —1, and the 
eigenvalues must have reciprocal moduli. In the sym- 
metric phase, both eigenvalues are unimodular, whereas 
the broken symmetry phase corresponds to the |Ai| > 
1,|A2| < 1 case. The exceptional point occurs when 
b — c — ±2, and there is a single eigenvector with eigen- 
value A = e^^^^^. Both the eigenvalues and amplitudes 
1^1 ^2 meet and bifurcate at the exceptional point, similar 
to the P7~-breaking transitions which occur in the eigen- 
value spectra of PT-symmetric Hamiltonians [THS] . 
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FIG. 5: Test of the criterion (171 for PT-symmetry breaking 
points of a three-layer hetero junction structure. The blue 
lines represent the eigenvalues of the S-matrix, which exhibit 
five symmetry-breaking points as the frequency lj is tuned 
over the selected range. The red line indicates the left hand 
side of Eq. (171. The heterostructure has a constant Re[n] = 
3, and the first and last layers are filled with gain and loss of 
Im[n] = ±0.005. The width of the central passive region is 
4% of the total length L. 

Each eigenvector of the S-matrix corresponds to a par- 
ticular choice of two coherent beams, simultaneously di- 
rected at each side of the heterostructure. The S-matrix 
transition can in principle be observed by tuning the 
complex input amplitudes, measuring the output ampli- 
tudes, and hence finding the scattering eigenvalues. One 
would actually need to do such "two-sided" interference 
experiments to detect the attenuating mode in the bro- 
ken symmetry phase, an interesting possibility which is 
currently being explored [23j . However, such experiments 
with two coherent input beams |17j are often inconvenient 
and difficult to perform. Therefore it would be preferable 
to have a criterion for the transition based on separate 
single-beam measurements. 

In Ref. [13] two such criteria were given for the phase 
boundaries in an arbitrary T^T-symmetric heterostruc- 
ture; however they both involve the relative phase of 
the reffection and transmission coefficients. One of these 
conditions is — = ±2it. Using the conservation 
relations ^ this can be shown to lead to the simpler 
condition |25j 



Rl + Rr 



-r = 1, 



(17) 



which involves only the transmittance and reflectances. 
The left hand side of Eq. ( 17 ) is greater than unity in the 



broken-symmetry phase and less than unity in the VT- 
symmetric phase. This provides a simple experimental 
criterion for locating the T'T-breaking transition point in 
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ID heterostructures. This criterion will be particularly 
useful if the quantity {Rl + Rr)/2 — T varies rapidly 
near the transition point. This appears to be the case for 
many heterostructures, as shown for example in Fig. [5] 
for a three layer structure. 
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FIG. 6: Reflectances and transmittance along the VT phase 
boundary for the ID PT-symmetric structure studied in 
Fig. [2|a), in the high frequency regime [ujL > 1725). The 
plots are given as a function of the gain/loss strength Im[n], 
while the frequency uj is simultaneously varied to maintain 
the system at the phase boundary. Vertical dotted lines indi- 
cate points where Rr = 0, for which Rl = 4 and T = 1 as 
predicted by Eq. ([T7|. 



Eq. 



(17) implies that for an ATR to coincide with 
the exceptional point, the non-zero reflectance must be 
exactly equal to 4, which is allowed but will not oc- 
cur without specific tuning. An example of such tun- 
ing is shown in Fig. |6] This plot is obtained by tuning 
both the gain/loss strength (Im[ri]) and the frequency, 
to keep the system along the phase boundary, and ob- 
serving the reflectances and transmittance. Two ATRs 
are found along the phase boundary. We note that a 
special set of solutions of Eqs. S and (17 1 are given by 
Rr, Rl,T = (p± l)^,p^, (p=F 1), where p is an arbitrary 
real number. Interestingly, the maxima of Rr,Rl,T 
in this simple geometry are given by this set of solu- 
tions with p = Re[n] in the high frequency regime where 
Im[7i] < 1. 



UNIQUENESS OF VT TRANSITION IN 
SCATTERING 



The generalized unitarity relations (|8])-(10) hold re- 
gardless of whether the eigenvalues and eigenvectors of 
the S-matrix are in the T^T-symmetric or T^T-broken 
phase; although the quantities in the generalized unitar- 
ity relations are related to the phase of the VT scattering 
system through the relation Eq. (17). There is, how- 
ever, some freedom of choice in the definition of the ID 
S-matrix, corresponding to permutation of the outgoing 
channels. The definition we used in Ref. [131 is given in 



Eq. ([g]), which is also widely used in mesoscopic physics 
[26] . In this section we will refer to the S-matrix defined 
in this way as Sq. In 5*0 the reflection coefficients are on 
the diagonal, and the outgoing channels are related to the 
corresponding incident channels by time reversal, which 
seems quite natural. In particular, the time-reversal op- 
eration T in this definition is represented by the complex 
conjugation operator. 

There is, however, an alternate definition: 



Sc 



TR t 



(18) 



which has also been used in the literature, including in 
one of the earliest works on T^T-symmetric scattering, 
by Cannata et al. [57]. This alternative definition of the 
S-matrix, which we will refer to as Sc, was subsequently 
used in the work on unidirectional invisibility of Lin et 
al. [in]. Because the permutation operation does not 
preserve the eigenvalues, these two different definitions 
of the S-matrix lead to different criteria for the sym- 
metric and broken symmetry phases, as well as for the 
phase boundary (exceptional points). This can lead to 
confusion, as well as raising questions as to whether the 
S-matrix eigenvalues and eigenvectors, and their transi- 
tions, are physically meaningful. 

Note first that both definitions lead to the same values 
for t, rji, tl, so they will give the same scattered state for 
the same input state. The issue is whether one or the 
other definition more closely reflects the phenomena of 
spontaneous VT symmetry breaking, as already known 
from Hamiltonian studies. In our earlier work on the VT 
transition in scattering systems |13j . we showed that the 
phase boundary of Sq corresponds closely to the anti- 
crossings of the poles of the S-matrix in the complex uj- 
plane (see also [15 ). The locations of these poles are 
independent of the definition of S; they reflect the in- 
ternal excitation frequencies of the scatterer, as well as 
the coupling of these excitations to the continuum. This 
suggested that the T^T-transition of Sq is indeed associ- 
ated with the VT transition of some underlying effective 
T^T-symmetric Hamiltonian. We have recently verified 
this point of view analytically and numerically, in collab- 
oration with others. The main part of that work will be 
presented elsewhere |2S] ; here we just state a few relevant 
results and show a numerical example corroborating this 
point of view. 

First, it is straightforward to show that the eigenvalues 
of Sc have the same general properties as those of 5*0, 
(even though they don't coincide). In particular, their 
product is —1 and they are either both unimodular or 
of reciprocal modulus. However the criterion for their 
exceptional points differs from that of Sq. Using a similar 
a,b,c parametrization of Sc as used earlier for Sq, one 
finds that the eigenvalues are given by: 



Al, A2 — 



1 



(19) 



Since both b and c are real, this expression shows that 



7 



when be > both eigenvalues are complex (and unimod- 
ular); whereas when be < 0, both eigenvalues are real 
and satisfy |Ai| — |A2|^^ ^ 1. Exceptional points oc- 
cur when b = or c = 0. From Eqs. ([7| and Q one 
sees that be = (l/T - 1) and so be > ^ T < 1 and 
6c < ^ T > 1, while b = 0(rL = 0) or c = 0{rR = 0) is 
the condition for T = 1. Thus each ATR is an exceptional 
point for Sc, and T > 1 corresponds to the "broken sym- 
metry" phase, whereas T < 1 to the "symmetric" phase. 
This is in contrast to 5*0 for which one has the criterion 
of Eq. (171 involving both T and the average of i?^ and 
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FIG. 7: Logarithm of the modulus of the eigenvalues of 5*0 
(red) and Sc (black) for the ID VT heterostructure studied 
in Fig. [2]^a); case (a) is with Im[n] — 3.995 x 10"'' and case 
(b) is with Im[n] = 4 x 10'^. |A| = 1 indicates the VT 
symmetric phase, and reciprocal values for |A| indicate the 
broken symmetry phase. Sc has multiple "transitions" spaced 
by the FSR and insensitive to Im[n]; 5*0 has a single transition 
which is highly sensitive to small changes in Im[n]. 

These two conditions for the transition and for the 
two phases of the S-matrix do not coincide (see Figs, [t] 
andjsj^a)) unless an ATR is tuned to occur at the phase 
boundary of 5*0 as we have shown in Fig. [6] We see that 
for this simple heterostructure, Sq has a single transi- 
tion to the broken symmetry phase (for a fixed Im[n]), 
while Sc has a series of transitions corresponding to en- 
tering and leaving the broken-symmetry phase in the 
high frequency regime (Fig. [t]). Each of these transi- 
tions begins at one of the two ATRs and ends at the 
other; thus the centers of the broken symmetry regions 
are spaced by the free spectral range of the unitary cav- 
ity. These "lozenges" of broken symmetry phase barely 
change when Im[n] is varied; Sc repeatedly enters and 
leaves the symmetric phase as we tune uj. In contrast, the 
single transition point of Sq moves substantially to lower 
frequency as Im[n] increases; once it enters the broken- 
symmetry phase, it never re-enters the symmetric phase 
at any higher frequency. This indicates that Sq, not S^ 
is measuring the breaking of VT symmetry. 



In Fig. [8] we show the decisive comparison. If we 
simply take the VT heterostructure shown in Fig. [2^, 
and impose Dirichlet or Neumann boundary conditions 
at the boundaries to the continuum, we have a non- 
Hermitian discrete eigenvalue problem with VT symme- 
try. Its energy spectrum (expressed as complex frequen- 
cies) makes transitions between real and complex conju- 
gate pairs (Fig.jsjb)), in a manner which perfectly follows 
the behavior of the eigenvalues of 5*0 and but not of Sc 
(Fig. [8](a)). Moreover, in Fig. [sjjc) we show the poles 
and zeros of the S-matrix; their symmetric distribution 
around the Im[A;] axis is a consequence of the VT sym- 
metry. Before the VT transition of Sq the poles have 
approximately the same value of Im[fc] as for the passive 
system, but just at the transition of Sq there is an anti- 
crossing in the complex plane and half begin moving to- 
ward the real axis and the other half recede further down 
in the complex plane [TS', 15^. For Im[fc] « 17 the system 
is very near the CPA-laser point for which a pole and 
zero coincide on the real axis. The eigenvalues of both 
Sq and Sc diverge/ vanish at this point because t^rj^ and 
all diverge at the lasing transition. Interestingly, for 
this value of Im[n], there are no ATRs after the lasing 
transition and T < 1 for all larger fc; the reasons for this 
are discussed in the Appendix. The same correspondence 
between the broken symmetry phase of Sc and the anal- 
ogous closed system hamiltonian holds for more complex 
VT hctcrostructures, such as that of Fig. ^c), where Sq 
has multiple broken phases . Thus we believe that at 
least for the ID case, there is a unique definition of the 
S-matrix, under which its VT transitions actually reflect 
the symmetry breaking in the underlying non-Hermitian 
Hamiltonian. 



VI. CONCLUSION 

We have derived generalized unitarity relations for the 
S-matrix of arbitrary ID T^T-symmetric photonic het- 
erostructures, including a conservation relation between 
the transmittance and the left and right reflectances. 
This conservation relation can be easily tested in experi- 
mental structures and used as a criterion of how precisely 
one has realized the VT symmetry. In addition, the con- 
servation relation leads to a simple criterion for identi- 
fying the exceptional point (s) at which the VT symme- 
try is spontaneously broken or restored. These excep- 
tional points are shown to be closely related to the VT- 
symmetry breaking transition of the underlying effective 
Hamiltonian of the system. 
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FIG. 8: (a) VT phase transition of 5*0 (red) and Sc (black) 
for the ID heterostructure studied in Fig. [2]^ a) but with 
n = 2 ± 0.208i. The broken symmetry phase of So is indi- 
cated by the shadowed area, (b) VT phase transition of the 
corresponding finite non-Hermitian cavity Ifamiltonian with 
Dirichlet boundary conditions. Its broken symmetry phase 
coincides with that of So. (c) Poles (red circles) and Zeros 
(black crosses) of the S-matrix (which are the same for So 
and Sc); anti-crossing of the poles occurs at the phase transi- 
tion of So. At kL ~ 17 the system is very near the CPA-Laser 
threshold point, where the cavity both emits laser radiation 
and perfectly absorbs the time-reverse of the lasing mode. 
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Appendix A: Properties of simple gain-loss 
hetero junction 



The VT heterostructure shown in Fig. [2][a), which 
consists of two uniform slabs of equal length and in- 
dex n, n* is the simplest example one can study of the 
class of ID VT symmetric photonic heterostructures, 
and it has been treated previously in [Hi [T^ . We will 
refer to this structure as the "simple hetero junction" 
(SH), and it is described by the transfer matrix Q with 
a = {a + a*) + i{P -h 7), 6 = -i{a - a*) + (7 -J), and 



a*) -I- (7 - where 

IcosAP n* 

a — 

2 2n 



sin AP, 



1 



2|n|2 
1 



[n* sin A cos A* -I- c.c] , 



7 = -[nsin A cos A* -f c.c], 



(Al) 
(A2) 

(A3) 



and A = nkL/2 is the complex optical path inside the 
left half. Since /3, 7 are real, so are b and c, and it is 
straightforward to check that ^ holds. This transfer 
matrix leads to certain simple properties. First, as men- 
tioned in the text, the SH has no trivial ATRs as we will 
show in subsection |A 1[ Second, below the VT symme- 
try breaking point it has many ATRs, roughly two per 
free spectral range of the passive resonator. Above the 
symmetry breaking transition it still has ATRs until it 
passes the lasing transition after which they disappear in 
the limit Im[nA:L] — )■ 00. We will discuss this behavior 
and contrast it with more complex heterostructures in 
subsection lA 21 



1. Absence of trivial ATRS 

The SH can be treated as having an air gap of vanish- 
ing width in between the gain and loss regions. Hence the 
absence of trivial ATRs is a consequence of the absence 
of reflectionless transmission resonances of such uniform 
amplifying or attenuating slabs in air. Below we first dis- 
cuss in general the transmission resonances of a uniform 
slab of refractive index n and length L/2 embedded in 
two semi-infinite media of index n; and Ur- 

For this simple setup the transfer matrix defined in 
Sec. II (( ^ ) = Af ( g ); see Fig. [T| takes the form 



1 



M = - 
2 



1 

1 - 



1 

nik 
riik 



cos A 



^ sin A 

infcsinA cos A 



1 1 

(A4) 

where n;, n^, n can be complex. A transmission reso- 
nance of an incident beam from the left side requires 
C = and 



A= l- 



cos A + i 



n 



sin A = 0. (A5) 



For the gain and loss regions in SH, when treated as being 
separated by an infinitely thin air gap, n; = = 1 while 
lm[n] = T 7^ 0. We immediately see that Eq. (A5) cannot 
be satisfied because sin A 7^ due to the finite imaginary 
part of n. This holds independent of n (r 7^ 0) and k{y^ 0) 
(see Fig. |9| . This finding is confirmed by calculating the 
reflectance directly (see Fig. [9]ja)). The same analysis 
can be extended to the slightly more complicated case 
shown in Fig. [2f^c), where all the ATRs are also found 
to be "non-trivial" as confirmed again by calculating the 
reflectances of the sub-units directly. We note, however, 
that trivial ATRs do exist in some VT structures. An 
example is the concatenation of even numbers of SHs. 
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FIG. 9: Reflectance of the gain (solid) and loss (dashed) re- 
gion in the heterostructure studied in Fig. [2][a). Two values 
of the refractive index are considered: n = 2 ± 0.2i (red), 
n = 2.5±0.1i (blue). Solid and dashed lines indicate Rl, Rr 
of the gain region and the loss region, respectively. 



2. Asymptotic properties of ATRs in the simple 
hetero junction 

In this main text we noted that ATRs for the SH 
disappear soon after the lasing threshold is passed in 
the broken symmetry phase of (see Fig. [lojja)). To 
understand this observation, we study the behaviors of 
Rfi and Rt in the large rkL limit. The reflection co- 
efficient from the "loss" side of the VT heterostruc- 
ture approaches the value given by the Fresnel formula, 
(1 — n)/(l -|- n), due to the suppression of interference 
effects by strong absorption. Surprisingly, the asymp- 
totic value of the reflection coefficient from the "gain" 
side, (1 -I- n*)/{l — n*), turns out to be the inverse of 
the Fresnel formula. This can be explained from the 
analysis shown by the inset of Fig. 10 'a). Due to the 



strong loss inside the "loss" side of the cavity, D' is 
given by the Fresnel relation C"(2n*)/(n* — n) = CV, 
implying \D\ — \exp[in*kL/2]D'\ = \r' exp[in* kL]C\ is 
much larger than |C| in the large rkL limit. Therefore, 
the scattering at the air-gain interface is approximately 
the time reversed process as if the "gain" side (which is 
the "loss" side in the time-reversed picture) were semi- 
infinite, i.e. with incident amplitude r* and reflected 
amplitude 1 in the air and transmitted amplitude D* , 
satisfying the Fresnel relation. The reflection coefficient 
in the original problem is then (1 + n*)/(l — n*). Thus 
RlRr — >■ 1 in this limit, which implies T — > from 
Eq. ^ and ATRs do not exist. 

In more complicated VT structures ATRs can take 
place in this limit. For example, the reflection coefficient 
connecting C" and D' approaches zero at the transmission 
resonances of the passive region in Fig. . The analysis 
above then breaks down and sharp changes of the trans- 
mittance and reflectances take place at these frequencies 
as shown in Fig. |10[ b). These transmission resonances 
through the passive region is a special set of solutions of 
Eq. (A5|. They require n; — n* and Im[n] = 0, and the 
transmission resonances occur at 



A = arctan 



2Im[n;]n 



(A6) 




FIG. 10: (a) Reflectances and Transmittance same as in 
Fig. ^sl) but at higher frequencies. Inset: Analysis of the 
reflection coefficient from the "gain" side. Letters with arrow 
represent the complex amplitude of the traveling waves at the 
nearest interface, (b) Reflectances and Transmittance same as 
in Fig.[2|c) but at higher frequencies. Inset: zoomed in on the 
two ATRs near ujL — 92.35. (c) Reflectances and transmit- 
tance of the central passive region in Fig. [2|c) placed between 
two semi-infinite regions of gain and loss with n — 2 ± 0.2i. 



The transmission resonances are given by Eq. ( A6 1 in which 
Im[ni] takes opposite signs depending on the propagation di- 
rection. Red curve represents T — 1 and its broken parts 
indicate T < 1. 



is small, these transmission resonances do not lead to 
ATRs of the whole heterostructure due to the multiple 
interferences taking place inside the gain and loss sub- 
units. In the large Im[n;]fc limit shown in Fig. |10| ^b), 
however, these multiple interferences are suppressed due 
to strong absorption/amplification, and ATRs arise from 
the resonances given by ( A6 ) . Note that these ATRs are 



Re[n;]2 — Im[n;]2 — n'^ 
In the frequency range shown in Fig. [2jc) where Im[n(]fc 



still "non-trivial" as the frequencies given by ( A6 1 are not 
the transmission resonances of the gain or loss sub-unit 
in the absence of the other. 

For the purpose of completeness, we mention a few 
more cases where transmission resonances of a single uni- 
form slab (i.e. the solutions of Eq. (A5)) exist. When 

, rir ,n can be treated as real (with negligible ab- 
sorption and no gain), two types of solutions of (A5) 
can be found. The first one is well-known, ni = rir, 
which requires sin A = 0; the second one is less well- 
known, n = y/niHr, which requires cos A = 0. It is 
easy to convince oneself that no other types of solu- 
tion exist for real indices. As one slowly increase the 
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gain or loss strength in the scattering layer, approx- 
imate transmission resonances can still be found, but 
their reflectances gradually increase and eventually be- 
come detectable. In Ref. [W a different case was studied 
where rir = l,Im[ri] — 0, Im[n/] 7^ 0. By noticing that 
cos A = cannot satisfy the above equation and tan A is 
real, Eq. ( A5 ) can be reduced to 



lm[niY 

tan A = —n 



{Re[ni] - l){n^ - 
Re[ni] - 1 



Re[ni]), (A7) 
(A8) 



It describes the transmission resonance from a loss/gain 
media to air through a passive slab, which gives rise to 
the novel "surface" lasing modes introduced in Ref. [25] , 
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